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Abstract: In order to achieve maximal stop mixing and m/j — 125 GeV, we consider extensions 
of minimal GMSB that include marginal MSSM-messenger superpotential interactions. Using a new 
approach to analytic continuation in superspace, we derive general formulas for the soft masses in the 
presence of such interactions, correctly taking into account the role of MSSM-messenger mixing in a 
general framework for the first time. We classify and catalog all possible such interactions consistent 
with perturbative SU{5) unification, and we survey the impact of turning on one interaction at a time, 
from the point of view of fine timing, spectrum and phenomenology. We find that the best models 
are fine-tuned to the sub-percent level and are accessible at the 14 TeV LHC. We highlight potential 
search strategies that can probe the characteristic spectra of these models. 



Contents 



1 Introduction 1 

2 A New Calculation of the Soft Spectrum 5 

2.1 Problems with the existing derivation 5 

2.2 A fresh approach to the calculation 8 

2.3 Formulas for type I models 9 

2.4 Formulas for type II models 10 

3 Models 11 

3.1 Type I Higgs couplings 14 

3.2 Type I squark couplings 16 

3.3 Type II couplings with mixing 18 

3.4 Type II couplings without mixing 20 

4 Phenomenology 22 

4.1 Type I squark models 23 

4.2 Type II models 24 

5 Conclusions 25 

A A Detailed Study of the QU^ Model 26 

A.l Applying our general formulas 26 

A. 2 Calculation using SUSY correlators 27 

A. 3 Analytic continuation method 28 

B Fine- Tuning Measure 30 



1 Introduction 

Recently, both the ATLAS and CMS experiments have announced the discovery of a Standard Model- 
like Higgs with TO^ « 125 GeV [1, 2]. This exciting result provides interesting hints and challenges 
for physics beyond the Standard Model. Supersymmetry (SUSY), long the preferred solution to the 
hierarchy problem, is highly constrained by this value of the Higgs mass. In particular, in the MSSM, 
large radiative corrections from stop/top loops are needed for nih = 125 GeV (see e.g. [3-7]). These 
contributions can arise either through extremely heavy, unmixed stops (Ms = ^m^yrn^ > 10 TeV), 
or through lighter stops with maximal mixing [8-10]: 

At ~ VgMs and Ms > 1 TeV (1.1) 

By transplanting the stops from 1 to 10 TeV, the theory grows two orders of magnitude more tuned. 
Since such a tuned model has little hope for ever being observed at the LHC, we will focus on generating 
the light, mixed stops in this work. 
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Large A-terms are essential for obtaining a heavy Higgs with hghter stops. This presents a special 
challenge for models of gauge mediated SUSY breaking (GMSB) (for a review and original references, 
see [11]), which are strongly motivated by the SUSY flavor problem, but do not produce A-terms at 
the messenger scale. Large A-terms can be generated through RG running driven by a heavy gluino 
[6], but this requires a very large messenger scale and, again, reduces both the naturalness and the 
likelihood of observing any superpartners at the LHC. 

In this paper, we will instead study models which directly generate large At at low messenger 
scales through marginal superpotential interactions between MSSM and messenger superfields. We 
will focus on fully calculable models of perturbative messengers coupled to SUSY-breaking spurious. 
There are then two types of marginal interaction terms: MSSM-messenger-messenger couplings, and 
MSSM-MSSM-messenger couplings. We will refer to these models as type I and type II models, 
respectively. It is useful to further divide the type I couplings into two distinct subclasses, those in 
which the MSSM superfield participating in the interaction is a Higgs, and those in which it is a 
squark. 

It is also useful to distinguish between couplings which give rise to mixing between the MSSM 
and messenger fields and couplings which do not. A prime example of a model with MSSM-messenger 
mixing is the W = XQsUs^ model, where $ is a messenger with the same quantum numbers as 
Other similar examples include W ~ XHuQs^ ($ mixing with U3) and W = XHuUs^ ($ mixing with 
Qs). Such models have been studied by many authors in the literature, including [12-20]. We only 
consider mixing in type II interactions, as the messengers in type I interactions can always be charged 
with a parity symmetry to forbid mixing. 

Fig. 1 displays this classification of models - into type I Higgs, type I squark, type II with mixing 
and type II without mixing. This classification will form the basis for the results presented in this 
paper. A primary goal of this work will be to describe models within these different categories and 
their general features with regards to fine tuning and phenomenology. 

The main theoretical challenge facing calculable models for large A-terms is something that was 
dubbed the Ajrr? problem in [21] (see also [22] for a recent discussion in a more general context). In 
direct analogy to the ^jl/ problem, models for A-terms tend to also generate one-loop soft mass- 



Type I: 

MSSM-messenger-messenger 



type I Higgs 



type I squark 



Type II: 

MSSM-MSSM-messenger 



type II 
with mixing 



type II 
w/out mixing 



Figure 1. A diagram illustrating our classification of models. 
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squareds. Such large soft masses would be disastrous for naturalness and/or electroweak symmetry 
breaking (EWSB). As shown in [21] (following [23]), spurion models which avoid the Ajrn^ problem 
must be of the minimal gauge mediation (MGM) type: 

Wmgm - K^*^ (1.2) 

in which all mass-scales in the messenger sector (consisting here of $ and $) originate from a single 
spurion X with {X) = Af + d^F . All of the models that we consider in this paper are assumed to have 
a messenger sector of this form. 

Although ^-terms are commonly viewed as trilinear soft terms, they actually arise as bilinear 
terms between the MSSM fields and their F-components: 

C D Ah^fI^H^ + Aq,fI^Q3 + Au.fIUs (1.3) 

After integrating out the F-components, this becomes AtHuQsUs with 

At = vMh^, + Aq, + Au^) (1.4) 

In order to generate At, it is mandatory that either H^, Q3, or U3 participate in the direct interactions 
with the messengers. In [21], only iJ„ couplings to messengers were considered, since this automatically 
preserves minimal flavor violation (MFV), one of the best features of gauge mediation. However, these 
models suffer from the residual "little A/rn^ problem" [21]: by integrating out Fh^ in (1.3), an 
irreducible contribution Smj^^ = +A'jj^ is generated. Since Ah^ ^ 2 TeV is required for maximal 
mixing and nih = 125 GeV, this results in an irreducible fine tuning at the ~ 10^ level in these models. 

Motivated by these considerations, we will broaden the scope of [21] and survey the complete class 
of spurion-messenger models for large At, including also non-MFV couplings involving Q3 and U3. 
As expected, these squark-type couplings result in far less fine-tuning than the Higgs-type couplings, 
since they do not generate a little A/mj^ problem. For simplicity, and to minimize flavor changing 
effects, we consider scenarios with only a single coupling introduced between MSSM and messenger 
fields. This structure may seem unrealistic; however, it is technically natural and facilitates focus 
on the "best-case scenarios." By restricting the models to complete SU (5) messenger multiplets and 
requiring perturbative gauge couplings up to the GUT scale, only a finite and manageable list of 
couplings is permitted. Under these assumptions, there are a total of 31 couplings. The complete list 
of these can be found in table 1. 

For computing the soft masses in these models, general formulas were derived in [12] using wave- 
function renormalization and the technique of "analytic continuation into superspace" [24]. These 
formulas are cast in terms of anomalous dimensions and beta functions of non-holomorphic couplings, 
and they are sufficient for all type I and most type II models. However, one runs into numerous compli- 
cations when attempting to apply them to type II models with mixing between MSSM and messenger 
superfields (such as the QU^ model described above). The problem stems from crucial assumption 
made in [12] that the wavefunctions are continuous through the messenger threshold. This proves 
to be incompatible with standard conventions for the beta functions and anomalous dimensions (e.g. 
those in [25]). Attempting to substitute standard beta functions and anomalous dimensions into the 
general formulas of [12], as many in the literature have done, leads to incorrect results.^ In fact, we 

^While one can correctly apply the formulas of [12] in the mixed type II case by using non-standard beta functions 
and anomalous dimensions (see [16]), this has only been performed in one specific model, and standardized formulas do 
not exist. 
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have found these complications to be so insidious that nearly every paper that studies models with 
MSSM-messenger mixing contains erroneous formulas. 

Faced with many incorrect results and an enormous amount of confusion in the literature, we will 
devote section 2 to describing a new approach to deriving the soft masses using wavefunction renormal- 
ization. The main idea of our new approach is to obtain the wavefunctions by direct integration in a 
manifestly holomorphic scheme, rather than first utilizing the beta functions for the non-holomorphic 
couplings, as done in [12]. This results in a conceptually cleaner approach, which circumvents the 
difficulties of [12] and leads to fully general, correct formulas which can be applied uniformly to all 
models, whether mixing is present or not. In appendix A, we provide many checks of our results for 
the QC/$ model described above, and illustrate the complications that arise in type II models with 
mixing, using this concrete example. 

With the correct and fully general formulas for the soft masses in hand, we will compute the soft 
masses at the messenger scale for each model in table 1, and then investigate the parameter space 
where m/i = 125 GeV. For such points in the viable parameter space, we will survey the phenomenol- 
ogy and fine tuning. In order to quantify tuning in these models, we utilize a tuning measure. Aft, 
that is based on the Giudice-Barbieri tuning measure [26], but with a slightly unconventional choice 
of underlying parameters. This choice is made to ensure that the tuning measure captures all sensitiv- 
ities successfully, does not introduce artificial tuning, and assigns comparable weight to uncorrelated 
contributions which cancel against one another. A detailed discussion and definition of the tuning 
measure is provided in section 3 and appendix B. 

Wc find that with the nih = 125 GcV constraint, the least-tuned spectra can be accessible at the 
14 TeV LHC, but are generally beyond the reach at 8 TeV. This suggests an intriguing possibility: 
that the failure to find superpartners so far at the LHC could actually be a consequence of nih = 125 
GeV, rather than a separate issue. 

Here is a summary of our results for each kind of model: 

• We confirm that the type I Higgs models (some of which have been studied before in [15, 17, 
21, 27]) are indeed fine-tuned at the level of ~ 10'* because of the little A/rnj^ problem. The 
spectra of the best points across the different models are rather similar. However, the prospects 
for observing these models at the LHC are rather pessimistic because all of the colored objects 
tend to be quite heavy. Some spectra do possess a light wino which could be produced at the 
LHC, but this is the exception rather than the norm. These spectra all exhibit slepton co-NLSPs 
with roughly 250-450 GeV masses - within the range that an ILC could discover. 

• On the other hand, in the type I squark models, there is no little A/rnj^ problem, so these 
models are considerably less fine tuned. To the best of our knowledge, these models have never 
been studied in detail before. All of these models possess points of relatively low tuning, with 
ApT ^ 10'^. However, as these models are not manifestly MFV, they can in principle confiict 
with fiavor physics constraints from precision experiment. For the purposes of this work, we 
adopt an agnostic stance toward flavor (e.g. we assume perfect alignment) and simply aim to 
show that the type I squark models are promising from a tuning point of view. Discussions of 
fiavor physics in these models will be deferred for future work [28] . The LHC phenomenology of 
these models is more promising. Often, there is an accessible stop (and sbottom in Q3 models) 
with slepton co-NLSPs, sometimes with a bino between the two states. Alternatively, some of 
the regions of lower tuning have gluinos and squarks light enough to be produced at 14 TeV. 
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• There are five type II models with mixing, and they have been studied in various guises in [12-20] . 
In three of these terms the interaction is top-Yukawa-hke {QU^, H JJ^ and HuQ^), and in two 
of them it is bottom- Yukawa- hke {QD^ and QH^^). The latter receive no appreciable influence 
from mixing, but they are more constrained by tachyons and/or EWSB, so they manifest with 
poor tuning. Conversely, the top- Yukawa- like terms are especially interesting as At is enhanced 
by receiving contributions from two sources (1.4). This enhancement is so effective that the 
mixed HuU^ model contains the points with the least tuning out of any model studied in this 
work. Though inferior to the best regions of all type I squark models, the H^Q^ and QU^ 
models have regions with lower tuning than the type I Higgs models and the rest of the type 
II models. The regions of least tuning in the QU4>h^ model are experimentally excluded by 
existing LHC searches because they have light O (1 TeV) gluinos and first-generation squarks. 
The other top- Yukawa-like models have spectra with gluinos and first generation squarks which 
will be accessible at 14 TeV running. Other production avenues, especially promising at an ILC, 
include very light slepton co-NLSPs and light Higgsinos entering at several hundred GeV. 

• Most of the type II models without mixing (see [17, 29-31] for some examples) are extremely 
tuned. These models tend to either suffer from the same issues as the type I Higgs models, or 
introduce additional tachyons which place the model in tiny corners of viable parameter space. 
However, the UsDs^d model has regions of slightly lower tuning. Additionally, these models 
depart from the slepton co-NLSP phenomenology of the type I models and will often have a bino 
NLSP. Gluinos accessible at 14 TeV could lead to exciting phenomenology. 

Our paper is outlined as follows: In section 2, we discuss the complications which arise when using 
analytic continuation in superspace to derive soft parameters in models with MSSM-messenger mixing. 
We then present a new, completely general framework for deriving messenger scale soft parameters in 
the presence of any MSSM-messenger interactions. In section 3, we catalog and survey the parameter 
space of the 31 possible couplings. Details are provided in a series of subsections for type I Higgs, 
type I squark, type II mixed and type II unmixed models. We discuss the phenomenology of the 
models with the least tuning in section 4. Section 5 contains a summary and discussion of future 
directions. We devote appendix A to validating the formulas of section 2 through multiple methods. 
The fine-tuning measure utilized throughout this work is detailed in appendix B. 

Note added: while this paper was in preparation, [32] appeared which overlaps partially with 
section 3 of our work. This paper also creates a catalog of the different models of MSSM-messenger 
interactions for large At. We note that their classification scheme differs from ours; our catalog contains 
several additional couplings absent from their study, namely Q4>Q(j)j), Uc/jqc/jh^, U(j)E4'Tji QQ4'~dj ^-s 
well as all couplings containing an SU{b) adjoint, ^24; and their formulas, being based on those of 
[12], neglect a proper treatment of MSSM-messenger mixing. 

2 A New Calculation of the Soft Spectrum 
2.1 Problems with the existing derivation 

In this section, we present a new method for calculating the soft spectrum induced by MSSM-messenger 
interactions via analytic continuation in superspace. The current state of the art (prior to this work) 
are the general formulas contained in [12]. Those formulas express the soft masses in terms of the 
beta functions and anomalous dimensions above and below the messenger scale. They are meant to 
be completely general, but are prone to misapplication whenever there is mixing between messenger 
and MSSM fields. 
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To understand the issues with the derivation of [12], we need to first review some aspects of 
wavefunction renormalization in the presence of operator mixing in SUSY theories. Let us define the 
theory above the messenger scale to be: 

K = z+m\-^, 

1 1 (2.1) 

Here t = log /i is the RG scale; since we are in the holomorphic scheme, the superpotential couplings 
do not run. The indices i, j, ... run over all messenger and MSSM fields (transforming in SM irreps). 
We allow for the possibility that any of the fields can mix with any of the others by taking the 
wavefunctions to be a general Hermitian matrix Z^j. 

In the holomorphic scheme, passing through the messenger scale is trivial, and the wavefunctions 
for the MSSM fields remain continuous. Below the messenger scale, the theory is of the same form, 
but without any messenger fields: 



b 



1 (2.2) 

W = -Xabc^a^b'^c 

6 

where now the a, 6, ... indices are summed only over MSSM fields only. Below the messenger scale, 
depends on M through its boundary conditions, so we denote this by Z~^^{t; M). 
The soft parameters are given by analytically continuing M — )• \/X*X and substituting {X) = 

M + e'^F into 

j d^eZ-,{t- \X\)<^i<^, D F {dxZ-,)Fl^<^, + h.c. + \F\\dxdx'Z-^)<^i<^, (2.3) 
so, to leading order: 



Aab - FdxZ^ 



b 



m 



2 _ |77|2 



- -,F\'dxdx^Z-, + \F\'dx-^Z-^dxZ-, 



(2.4) 



with the derivatives evaluated at i = log \X\. The second term in m^j, comes from integrating out the 
F-components of Note that because these expressions exist below the messenger scale, all indices 
correspond to MSSM fields only. 

Finally, the derivatives of the wavefunctions can be obtained from the integral expressions: 

Z-bii-^M)= dt''^+ dt'"^ (2.5) 

JlogAf JlogAuv 

and using the relation between the wavefunctions and the anomalous dimensions:^ 

f = -2VhV (2.6) 



^We note that [12] also used an incorrect equation for wavefunction renormalization, ^ = 7Z. From this starting 
point, [12] derive a formula for that contains a spurious extra term which goes as the commutator of 7 above 

and below the messenger scale. This must be incorrect, because it is not Hermitian. In any event, it drops out of the 
Lagrangian when sandwiched between ^a^b- 
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where Z — V^V, and 7 is the matrix of anomalous dimensions as a function of the non-holomorphic 
running couphngs, which we wih denote by X{t). These are related to the holomorphic couplings via 



kjkit) = X^'rk'{V'')^'^iV~')rJiV^')k'k (2.7) 

One can check that with these equations, the beta function for Xijk{t) is the standard one given in 
e.g. [25]. 

The problem in the derivation of [12] arises at this stage. Chacko & Ponton view (2.6) as an 
equation which determines ^ in terms of the non-holomorphic running couplings Xijk{t). These are 
defined in the non-holomorphic basis where the Kahler potential is canonical. If these were the only 
running couplings, then the derivation of [12] could be safely applied. However, this is not the case: 
in addition to Xijk running, the non-holomorphic couplings kij also run. That is, analogous to (2.7), 
we also have: 

~n,,it)^K.,,,{V-')MV-%, (2.8) 
Thus, the theory in the non-holomorphic basis, where the derivation of [12] takes place, is actually: 

K = $11., 

1. 1 . - (2.9) 

W = -A,jfc(i)$,$j$fc + -K„(t)X$,$j 

Even with (t) entirely aligned with the messenger directions in the far UV, the presence of MSSM- 
messenger mixing will cause it to become nonzero in the MSSM direction by the time we reach the 
messenger threshold. Now, integrating out the messengers will introduce some dependence on the 
MSSM fields, schematically: 

^^ aa^a (2.10) 

Substituting this back into (2.9), we obtain the theory below the messenger scale where the effective 
Kahler potential for the MSSM fields shifts discontinuously, 

K^K + ^^C^C^ (2.11) 

This potential is neither canonical nor continuous through the messenger threshold. 

To summarize, the problem with applying the derivation of [12] in the presence of MSSM-messenger 
mixing is that it relies on the non-holomorphic scheme where the Kahler potential is canonical and the 
superpotential couplings run. However, it also assumes that the wavefunctions of the MSSM fields are 
continuous through the messenger threshold. Together, these assumptions prove incompatible with 
the standard beta functions for the non-holomorphic couplings. 

The discontinuity in the wavefunctions (2.11) amounts to an additional contribution to the soft 
masses in the presence of MSSM-messenger mixing, which is missed in the formulas of [12]. One could 
attempt to take this into account by fleshing out the line of reasoning above, as done in appendix A 
for a particle model. Alternatively, one can perform an extra unitary rotation to undo the effect of 
the RG (2.8) and to eliminate the extra contribution (2.11). While this latter method yields correct 
results (it was done for a specific model in [16]), a general formula derived from this procedure does 
not exist. Furthermore, the price one pays in this approach is that the anomalous dimensions and beta 
functions are no longer given by standard formulas; in particular the matrix of anomalous dimensions 
is no longer Hermitian. 
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2.2 A fresh approach to the calculation 

In this paper, we will take a fresh approach to the problem of computing wavefunction renormalization, 
one that will overcome the problems discussed above. As reviewed in the previous subsection, in 
the standard implementations of analytic continuation, one first puts the Kahler potential into a 
canonical form, solves the beta function equations for the non-holomorphic couplings, substitutes 
these into formulas for anomalous dimensions, and integrates these to get the wavefunctions Z which 
were canonicalized in the first step. This methodology seems rather convoluted, since the running 
of the non-holomorphic couplings is nothing other than wavefunction renormalization. Wouldn't it 
be conceptually simpler to stay in the holomorphic basis and directly integrate a single differential 
equation for Z, rather than integrating two differential equations which are really expressions of the 
same underlying physics? 

All that is required is to view (2.6) as an equation for Z itself, and not as an equation which 
determines Z in terms of the non-holomorphic couplings. At one-loop, the anomalous dimensions are 
given by: 



1 n,kr. 



^ l&TT^^"'^^ {2'^'^'^ ^i'l'^^km ^en* " 2cl. Z^> j, g'^j Vj,j 



(2.12) 



Here df^ is a standard multiplicity factor present in the one-loop anomalous dimensions; roughly 
speaking it counts the number of fields of type k and £ that can talk to field i through the interactions. 
Concrete examples of df^s will be given in later sections. In the second line of (2.12), we used the 
facts that d^^ is a function of only the gauge representations of i, k and £, and that Zij and Vij can 
only connect two fields in the same representation. Substituting (2.12) into (2.6), this becomes 

^ = [Z{t)-\ g{t)] ^ (^rff A:,,Z-,i*Z7„i* A„„„ - 24Z,,,g2^ (2.13) 

The y's have disappeared, and we have the desired form of the differential equation for Z given in 
terms of Z itself. 

It remains to compute the first and second derivatives of Z~f^ with respect to X and X* . Here we 
can follow essentially the same steps as in [12]. Repeatedly using (2.5) and (2.13), we arrive at: 



5Z-,(logAi;|A| 



dX 

32 ' 



2A 



dXdX* 



4|A|2 \^ \dZ,J dZr^ V 5.g. ; ^'^ dgV 



(2.14) 



where A(. . . ) means the discontinuity of (...) across t = log |A|. These results are clearly analogous 
to the formulas in [12], but they are more broadly applicable. In particular, there is no complication 
to applying this formula to models with MSSM-messenger mixing. Substituting the explicit formula 
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for Gij given in (2.13) into (2.14) and then into (2.4), we obtain at leading loop order, 

<5™a6 = 255^4 iKik^bjk) {KerrA*jhn)^ - (Kik^bjk) A {Xi£m>^jlrn)^ (2.15) 

A (a:,, -A^,) A (A:,,A,fe,) - dl^Crg'rA (A^ -A,.,)) 

where A = F/M, and C;^'= = cj. + 4 + is the sum of the quadratic Casimirs of each field interacting 
through Xijk- In this expression, we have not bothered to write the usual GMSB term (hence the 5 
in front of m^;,), which comes from the last two terms in the second line of (2.14). All indices are 
summed over except for a and b. As discussed below (2.4), the sum over c in the second to last term 
must include only MSSM fields. 

This is our final, general result for the one-loop A-terms and two-loop mass-squared terms induced 
by MSSM-messenger interactions. On top of the standard GMSB contribution to the soft masses, one 
must add to this expression an additional term appearing at one-loop, but suppressed by ^ [21]: 

Sml,^i-ioop - -^^^^<^ A (A^ .A,,,) ^ (2.16) 

where 

h{x) = ^{{x- 2) ln(l -x)-{x + 2) ln(l + a;)) = 1 + + O (x") (2.17) 

Using these formulas, one can derive correct expressions for the soft masses and A-terms for any model 
with any number of type I and type II MSSM-messenger interactions, including those with mixing 
between any and all sectors. However, in the following subsections, we present specific simplified 
formulas for models containing type I couplings only and type II couplings only. These formulas are 
used throughout this work. 

2.3 Formulas for type I models 

In the previous section, we derived the most general formulas for the MSSM soft masses, including 
the possibility of arbitrary mixing between MSSM and/or messenger fields. Now, we would like to 
specialize to models which involve only type I or type II couplings, beginning with the type I case. 

To improve the readability of the formulas, it will be convenient to introduce messenger-only 
indices A, B, . . . . The indices a, 6, ... will continue to run over MSSM fields only; and i, j, ... 
will continue to run over all fields. Finally we will denote MSSM-messenger interactions with A, but 
MSSM-only interactions (the usual Yukawa couplings) with y. 

In the type I models, the interaction is of the MSSM-messenger-messenger type: 

W = ^XaBC^a^B^C (2.18) 

In these models, one can always impose messenger parity, so that MSSM-messenger mixing does not 
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occur. Specializing (2.15) to the type I case, we obtain: 

A-ab = -l^TT^d^^ KiBC^bBC^ 
'^'^afc = 2567r* iy^^ KiBC^bCEKBoKoE + '^d^'^' d^^ Xlg^XcBcKoE^bDE (2.19) 

" l^d^^d^^ yl^^ybde^cBcKBC " ^r^^ glK.BC'^bBC^ 

If we further speciaUze to the case of no MSSM-MSSM mixing and no messenger-messenger mixing 
(e.g. only a single coupling between MSSM and messenger sectors), then this becomes 

(^"^a = 2567r4 {^^^ '^'B'\^aBc\^\'XcBD? + ^ rff^^f ^ | AaBC H ^aDB P (2.20) 

which agrees exactly with the formulas given in appendix A of [21]. 
2.4 Formulas for type II models 

In our framework, it is clear that the type II models with mixing are really no different than the type 
II models without mixing. When a particular model does not have any MSSM-messenger mixing, 
some of the terms in Smf^ are simply zero. Nothing special needs to be done and the mixing is fully 
accounted for by the formulas without requiring any further treatment. In all type II models, the 
interaction is of the MSSM-MSSM-messenger type: 

W = ^Aa6C*a$6*C (2.21) 

Specializing to this case, the soft SUSY breaking terms are now given by: 

^ab = ^Yg^^a^^'^C'SAbcsA 

^l^lb = 255^4 ( ^*acB^bcC>^deB>^*deC + d'f df^ Xl^^^XbeBKdcKleC 

cB * cd B . * 1 cS e * * ^^'^^^ 

+ da^ df^ X^^gXceBKleC^bdC " d"^ yl^^VbdeKf BKf B + Ij^d"^^ d^J ycefV*defKicB^bdB 

+ 2'^''adl^y*acdycef'XbdBKfB + l^d'f d""^ yl^dVdef^bcBKf B " 2(iJJ^C°''^C/^A*cB AfccS^ A^ 

The first two terms in the last line are the additional contributions which arise from couplings with 
MSSM-messenger mixing. It is easy to see that these indeed vanish if there is no MSSM-messenger 
mixing, since in that case XefB and ycef cannot be simultaneously nonzero. It can be checked that 
naively substituting the standard beta functions and anomalous dimensions into the formulas of [12] 
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misses precisely these extra terms. 

If we assume there is no MSSM-MSSM or messenger-messenger mixing, then (2.22) becomes: 

^ ^ 1 (2.23) 

+ dfdf'\XacB?\XadC? " «^ lyacrfPl A./B P + l^df \y,,f\''\XacB? 

+ \d''^d!'J VlcdVcefXadBXlfB + ^d^^ d^J y^^^ydef X*acB>'ef B " 2d^^C;'=-^g^ | AacB 

where, as before, the first two terms in the last line will vanish in the absence of MSSM-messenger 
mixing. 

3 Models 

From the general considerations of the previous section, we now turn our focus to surveying the 
different types of MSSM-messenger interactions. To produce our catalog of models, we impose a few 
conditions. First, we require that the messengers come in complete, vector-like SU(5) representations 
and that the SM gauge couplings remain perturbative up to the GUT scale. The relevant SU{5) 
representations and their SU{3) x SU{2) x U{1) decompositions are: 

<^5 ® 05 (3, 1)_ 1 ® (1, 2) 1 e (3, 1) 1 (1, 2)_ 1 

' - ' _ (3 1) 

(1)10 © -^10 ^ (3, 2) 1 © (3, © (1, l)i © (3, 2)_ 1 © (3, 1) 2 © (1, 

(^24 ^ (8, 1)0 © (1, 3)0 © (1, 1)0 © (3, 2)_ 5 © (3, 2) 5 

6 6 

In order to keep the gauge couplings perturbative up to the GUT scale, the total messenger contribution 
to the beta function must satisfy |A6| < 6. 0i, ^5 © cAg-, (^10 © (l>To and (j)2A contribute 0, -1, -3 and -5 
to the SU{5) beta function, respectively. 

Second, we will only consider models where a single superpotential coupling A is turned on between 
the MSSM and messenger fields. For type I models, our interaction superpotential will consist of a 
single term: 

W = XqY,<l>^'|)^ (3.2) 

1=1 

Here, the MSSM field, q, must be one of either Hu, Q3, or U3 in order to produce large At- Meanwhile, 
(\)\ are messenger fields transforming in S'C/(3) x SU{2) x U{1) irreducible representations. In (3.2), 
we have included the possibility that the messenger number, Nj^, may differ from one - a viable option 
for type I models. Note that we have taken all couplings to the multiple pairs of messengers to be A 
for simplicity. 

For type II models, the story is nearly identical. Here, our interaction superpotential will be: 

W = Xqq'cj) (3.3) 



- 11 - 



with q, q' being MSSM fields (at least one of which must be either Q3, or U3), and (j) a messenger 
transforming as an SM irrep. Note that for type II models there is no option to talk to multiple 
messengers; the messenger number is always one for the contribution from the MSSM-messenger 
interactions. The messenger number for the GMSB contribution could be greater than one, but this 
always serves to make the model more fine-tuned (it will reduce At /Ms), so we will restrict ourselves 
to Nm = 1 for type II models. 

Under these constraints, there are 31 models in all - 15 of type I and 16 of type II. These are 
cataloged in Table 1. Each model in this table can be parametrized in a uniform way. As in [21], we 
will choose the parameter space to be: 

X, ^, N^, a) (3.4) 



M' 

For a given value of the first three parameters, increasing A = F/M simply increases the overall 
scale of the sparticle masses and A-terms (while keeping the stop mixing fixed). Since to/j increases 
monotonically with A (by having heavier stops running in the loops), A is uniquely determined by 
imposing nih = 125 GeV. This procedure can fail if some sparticles are tachyonic or if the basic 
conditions for electroweak symmetry breaking cannot be satisfied. 

In almost all of the models which we discuss, either U3 or does not talk to the messenger 
directly. That field will run tachyonic at high values of A due to the leading terms in the soft masses 
being of the form, 

{l6n'Y'^ = ~B,y^X' + B,gt. (3.5) 

where Bi, B2 > are the contributions from the MSSM-messenger coupling and the standard GMSB 
contribution, respectively. So at a large enough value of A, the squark becomes tachyonic. Solving for 
A, this happens near, 

A™..^^/|5|«yj (3.6) 

This squark tachyon ceiling appears in nearly all models, although RG running from the messenger 
scale perturbs the actual value of Xmax (which manifests in gradual curves, as opposed to horizontal 
lines in (A, ^) space). 

To explore each of these models, a dense grid of points in A and ^ is generated, with A chosen to 
satisfy the nih = 125 GeV condition. We use SOFTSUSY [33] to run the spectra from the messenger 
scale to the TeV scale. As in [21], contour plots in the (A, plane completely characterize the viable 
parameter space of the model. If a particular point gives a valid solution, then the fine-tuning at that 
point is calculated using our tuning measure. As discussed in more detail in appendix B, the tuning 
measure is defined as 

d log m'^ 

Aft = maxAi where ee — (3.7) 

d\ogAf 

where A^ G {Aa, A3, Af , Ai_;oop, m} with Ax = A^A, A3 = g^A, At = yjA and Ai^ioop measures the 
variation of the one-loop term (2.16). Each quantity is treated independently and varied separately. 

The least tuned point located in each model is cataloged in table 1. The spectra of the models 
with low tuning will be discussed in detail in section 4. The UHu4>q type II model is the least tuned 
{ApT ~ 850) out of all models, however, the rest of the type II models fare poorly with regards to 
tuning. In general, all of the type I squark models enter with a relatively low tuning measure of 
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# 


Coupling 


\Ab\ 


Best Point \^,X} 


\At\/Ms 


Mg 


Ms 


ImI 


Tuning 


I.l 


Huft'^ r (t>l S 




{0.375, 1.075} 


1.98 


3222 


1842 


777 


3400 


1.2 






{0.25,1.075} 


1.99 


3178 


1828 


789 


2450 


1.3 




4 


{0.25,1.3} 


2.05 


2899 


1709 


668 


3200 


1.4 




4 


{0.125,0.95} 


0.58 


11134 


8993 


2264 


4050 


1.5 


f/'^ /- f/'24 S 


6 


{0.225,1.000} 


0.54 


13290 


9785 


3408 


3850 


1.6 


O?^ r (Z>24 W 


6 


{0.15,1.025} 


0.67 


11835 


8637 


3259 


3410 


1.7 


H^ihw ri4'2i X 


6 


{0.3, 1.425} 


2.04 


3020 


1743 


576 


3500 


1.8 


0<t>Tn 7=;<i'l 9 


3iV„ 


{0.534,1.5} 


2.82 


4336 


1274 


2056 


1015 


1.9 


Q'h D'h L 




{0.353,0.858} 


2.67 


4247 


1342 


2058 


1015 


I.IO 


Q4'W,u4>5,H^ 


4 


{0.51,1.788} 


2.65 


4040 


1318 


2301 


1275 


I.ll 




4 


{0.378, 1.245} 


2.76 


4020 


1257 


2292 


1260 


1.12 


Ud>Tn Tf4>^ s 




{0.476, 1.622} 


2.62 


3815 


1347 


2070 


1030 


1.13 






{0.301,0.908} 


2.91 


3829 


1199 


2061 


1020 


1.14 




4 


{0.37,1.352} 


2.81 


3575 


1220 


2312 


1285 


1.15 


U (biQ e4>'^ ~n 


4 


{0.51,1.972} 


2.63 


3526 


1312 


2310 


1280 


II. 1 




1 


■To 55 1 641- 


2.02 


769 


1965 


2738 


1800 


II. 2 




3 


10 009 1 0671 


2.14 


2203 


1628 


543 


850 


II.3 


CJlli, d)'] n TT 


3 


{0.269,1.05} 


2.27 


2514 


1458 


439 


1500 


II.4 




1 


{0.37,1.2} 


1.78 


2597 


1829 


3553 


3020 


II.5 




1 


{0.15, 1.19} 


1.45 


2497 


2108 


3773 


6050 


II. 6 




1 


10 45 ll 


0.22 


7943 


9870 


3610 


5000 


TT 7 


TTDAy- 


1 
1 


if) 91 1 9f\\ 


9 m 


lO 1 ^ 


1 ^'?4 

lOO^ 


9QQ8 


91 ^0 


TT S 
11. o 




1 
1 


rn 1 A 1 9T. 


1.01 


lOUl 


1 904 
1ZU4 


zzuo 


'i7nn 

<J t uu 


TT O 

ii.y 






|U.440, i.4Dj- 


i.oy 


zUU4 


1 1 ou 


OO ( o 


z ( oU 


11.10 




tj 




2.13 


9943 


1 64Q 


282 


3500 


11.11 


HuL4'l,S 


1* 


{0.15,0.675} 


0.54 


7103 


8166 


3714 


4930 


11.12 




5 


{0.296,0.96} 


0.53 


12629 


9660 


3333 


3780 


11.13 


HuL(t)24,W 


5 


{0.212,0.96} 


0.65 


11487 


8710 


3687 


3380 


11.14 


HuHd<Pi,s 


1* 


{0.125,0.675} 


0.55 


7049 


8051 


3255 


5000 


11.15 


HuHd4>2i,s 


5 


{0.20, 1.00} 


0.57 


12047 


9213 


1628 


4220 


11.16 


HuHd(j)24,W 


5 


{0.2,0.946} 


0.64 


11571 


8789 


3665 


3460 



Table 1. All possible marginal MSSM-messenger couplings compatible with a perturbative 5(7(5) framework 
are tabulated here. The point with the least tuning in each model is also presented. The tuning measure 
used is defined in (3.7) and is discussed more in Appendix B. Additionally, the values of \At \ /Ms, Mg, Ms 
and at this least tuned point are shown. Models with \At\ /Ms < 1 rely on heavy stops as opposed to 
mixed stops. Models 11.11-13 generate large neutrino masses. Models 11.14-16 possess a ii/B^ problem. In 
the third column, |A6| refers to the messenger contribution to the SU{5) beta function. As the singlet does 
not contribute to GMSB, models 11.11 and 11.14 are assigned an additional (j>s ffi (I'd- 

ApT ^ 10"^. Many of the models involving Higgs fields have very large Ms (and small \At\/Ms) 
because they are relying on heavy stops to generate — 125, as opposed to using maximal mixing. 
As these models are unable to achieve maximal mixing without substantial tuning entering elsewhere 
(due to the little A/mj^^ problem), we make no effort to optimize the tuning in these models by 
scanning regions of parameter space where the MSSM-messenger contributions are small. Details 
concerning the various models will be discussed in the next subsections. 
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Aft for Hu(/)Ha(/>s Apx for HuCpQCpu 




A/M A/M 



Figure 2. Examples of the type I Higgs model parameter space. Left: HuXp^ jj^(jii,s with Nm = 6. Right: 
Hu4>io,q4'io,u with Nm ~ 2. In these, and all other type I Higgs models, several features emerge. At high 
slepton tachyons appear. At low jj, increasing A causes issues with EWSB. Moderate jj allows for the largest 
A, but this growth is truncated by a tachyonic stop, the scale of which is defined in (3.6). At moderate A and 
low i.e. the upper left corner of the viable region, both plots display where most type I Higgs models have 
their region of least tuning. 



3.1 Type I Higgs couplings 

We now survey the models in some detail, beginning with the type I Higgs models. As discussed in 
[21], these models are all MFV, but they have high tuning because of the little A/mjj problem. 

With only a single MSSM- messenger coupling of the form AiJ„(/)i(/)2 appearing in the superpoten- 
tial, the general form of the soft parameters from (2.20) can be easily specialized to the case of type I 
Higgs models: 

5m% = ^—rdnX'^yfA'^ 

^ 2567r4 " ^* 

2 

^ 2567r4 ^* 

where the multiplicity factors are defined dn = df^^^ , = d^"'^^ + d^^"^^ and Cr = c^" + cf^ + cp 
is the sum of the quadratic Casimirs. The values in each model for dj^, d^, and Cr are displayed in 
table 2. 

Qualitatively, all of the type I Higgs models have nearly congruent parameter spaces, with three 



487r2 



M2 



(3.8) 



- 14 - 



# 


Model 




I.l 

1.2 
1.3 
1.4 
1.5 
1.6 
1.7 


Hu4'w,q4'w,u 


3 (^,§,0) 

SAT (13 3 8) 

3 3 (J- 3 8\ 

V 30 ' 2 ' 3 / 

1 3 (^,2,0) 

1 3 (^,§,0) 

2 2 V 10 ' 2 ' 

V 30 ' 2 ' 3 / 



Table 2. The definitions of the parameters appearing in the formulas for the type I Hu models. The 
multiplicity factors are defined iIh = dti^^ j "^0 = '^^"'^^ + '^'^^'^^ ™d = c^" + cf^ + c^^ is the sum of the 
quadratic Casimirs. 

common, characteristic features emerging, see fig. 2. These features were discussed at length in [21]; 
let us briefly review them here. At lower increasing A gives a very large positive contribution to 
77i|^ which quickly causes issues with EWSB. On the other hand, raising increases the negative 
contribution to from the one-loop term (2.16), which makes the viable region of A larger. These 
two features combine to form the left, slanted edge of the viable regions shown in fig. 2. Increasing -jj 
even further eventually caps A by confronting it with a stop squark tachyon and drives even more 
negative (by raising the 1-loop term). As the model approaches ~ 1, eventually the increasingly 
negative m^j drives the sleptons tachyonic before the EWSB scale. 

In different regions of the parameter space, the tuning parameter is dominated by different quan- 
tities. For low values of A, the tuning is set by /i and A3. This is the case in all models, because 
as A — > 0, the contributions from both A^ and Ai_/oop vanish and At is generally subdominant. For 
moderate A and small -p^, A^ is the biggest contributor to tuning. It is here that we find the lowest 
values of /i, and hence the lowest values of the tuning parameter. Finally, for larger the Ai_;oop 
contribution governs the tuning. 

Finally, let us comment on some of the differences between type I Higgs models that are apparent 
from table 1. 

• The first two models - Hu4^ h^4'i,S a-nd -ffM0io,Q0io,c/ ~ were studied in detail in [21]. For these 
models, different are possible. In type I Higgs models, increasing the messenger number 
within a specific model decreases the tuning for that model. This happens because At and m? 
both scale as the number of messengers, thus At oc ^NmMs, which makes it easier to achieve 
At ^ \/6Ms for a maximal contribution to the Higgs mass. For all the other type I Higgs models, 
only Nm = 1 is possible because anything greater would violate the |A6| < 6 constraint. 

• Models 1.4, 1.5 and 1.6 have small At/Ms- These models forfeit maximal mixing in exchange for 
heavy stops. Models 1.4 and 1.5 are both identical to model I.l in terms of their A-terms and 
their A contributions to soft masses, but in terms of their GMSB contribution to soft masses they 
have effective messenger number 4 and 6 respectively. As these models are unable to achieve 
maximal mixing without substantial tuning entering elsewhere, we do not optimize the tuning 
further by scanning into regions of parameter space where the MSSM-messenger contributions 
are vanishing. 
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• Models 1.3 and 1.7 both have shghtly better tuning than 1.4-6. This is because these models 
receive an enhancement from the multiplicity factor dn^ — 3. This enhancement provides larger 
A-terms and allow significant stop mixing to be achieved. 

3.2 Type I squark couplings 

In the type I Higgs models, received a large correction +74|^ , leading to the little A/toI^ problem 
and greater fine tuning. Type I squark models receive an analogous correction to mg^ or rnf^^, but 
this poses much less of a problem, as the electroweak symmetry breaking scale is only sensitive to 
the stop masses at loop level. Without the little Ajrr?^^ problem, the type I squark models exhibit a 
significantly reduced tuning with respect to the type I Higgs models. 

While the type I squark models fare well with regards to tuning, the lack of an MFV structure 
makes them more dangerous with regards to flavor. However, the contributions to flavor changing 
processes only arise from A-terms, which are slightly safer than the off-diagonal sfermion mass terms. 
Some of the most dangerous flavor changing process involving A-terms come from EDMs and q q'^ 
processes. However, if there is a very strong alignment between the MSSM Yukawas and the ^-terms, 
then there are never any flavor physics issues arising from the A-terms simply because the coupling does 
not change flavor. Obviously, for a sufficiently small perturbation around this limit, flavor constraints 
will still not be violated. Precisely how small this perturbation must be and whether this alignment 
can be achieved naturally are interesting questions for future studies [28]. Regardless, while flavor 
appears alarming in these models, these concerns are insufficient to invalidate the models outright. 

With only a single MSSM-messenger coupling of the form XQ(f>i(j)2 or XU(t>i4)2 appearing in the 
superpotential, the general form of the soft parameters can be derived from (2.20). For the type I 
Q-type models, we have. 



Ao = ^doA^A 

W 167r2 



(3.9) 



2567r4 ^ ^* "-^ 2567r4 

2 2 

(5mr; = jdnX^yih? Sml, = jdnX^ylly^ 

^ 2567r4 ^* ^ 2567r4 



and for the type I J7-type models, we have. 



Au - --^duX^A 



"^"^^^ = -25^^^^^'^*^' = -25^^'^^^'^*^' 

where as in the Higgs type models, we define the multiplicity factors dg = dg'^^ and d(j = df/-'^^ , 

d(j> = d^'j'^ +(^'^2^ or = d.'^i^ +'^02'^^ Cr — c^'^ +cp +cp as the sum of the quadratic Casimirs. 
The values of dq^u, d^, and C'r in each model are shown in table 3. 

As was the case in the type I Higgs models, increasing Nm improves the tuning, however, here 

1 

At (X NmMs, due to one stop scaling as m~ cx N'^ and the other as to? oc Nm. The Nm enhanced 
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# 


Model 


dg dcf, Cr 


# 


Model 


djj d^ Cr 


1.8 
1.9 
1.10 
1.11 




( l30' 2' 3/ 

JV„i o V3O' 2' 3/ 

1 fis 3 8^ 
^ V3O' 2' 3/ 

2 6 (^,§,4) 


1.12 
1.13 
1.14 
1.15 


U(t)w,g4'5,H^ 


4 (A,o,|) 
2A^„, 4 (1,0,4) 
2 4 fis 3 8\ 

^ ^ V 30 ' 2 ' 3 / 

1 4 (1|,0,|) 



Table 3. The definitions of the parameters appearing in the formulas for the type I squark models. Q models 
are on the left, U models on the right. dQfjj) = d'g^^^ is the multiplicity factor, = d'^^^'^^ + d'^'^'^'^^ and 

Cr = c?'^' + 4^ + cp for the type I Q{U) models. 



models, 1.8-9 and 1.12-13, possess slightly lower tuning than the models which cannot capitalize on 
this feature. 

The parameter space of the type I squark models possesses a common, characteristic shape with 
two distinctive features - the "horn" and the "throat" - as shown in fig. 3. At high A, the squark 
which does not interact with the messengers (either U3 or Q3) is tachyonic, as discussed around (3.5). 
As ^ decreases, the messenger scale increases, so the stops have more time to run negative, thus 
forbidding smaller values of A and resulting in a curved viable region. Meanwhile, at intermediate A 
and larger -p., the squark which does interact with the messengers can go tachyonic, as its soft mass 
is given schematically by, 

where the Bi are positive numbers, and the function h{x) is defined in (2.17). The combination of stop 
tachyons above and below in A create a horn-like feature that extends into higher A and Going 
even lower in A, eventually (3.11) will rise again, and the stop will cease to be tachyonic. This results 
in an excluded interval at intermediate A, which we will call the "throat" region. It exists as long as 
^ is greater than some critical value: 



(3.12) 

MJ,r^t V 167r2B2S4 ^ ' 

where we have approximated h{A/M) w 1. 

At small A, the tuning is set by /i and A3. The tuning in the interesting regions of these models 
is everywhere a balance of /i and A^, and at higher values of A (i.e. on the horn region) the latter 
is leading. The region of lowest tuning in these models sits roughly on the underside of the base of 
the horn. This is sensible because one of the stops will run tachyonic near here, generating a larger 
At /Ms- At the base, the tuning is set by fi. For two models (1.9 and 1.13 - both shown in fig. 3), a 
second region of comparably good tuning sits in the middle of the horn. Although the stops are much 
heavier here, the tuning does not suffer greatly because the messenger scale has significantly decreased, 
and with it, \m'jj | has decreased resulting in a lower value of fj.. Other type I squark models exhibit 
some regions of parameter space with a similarly decreasing A, but the tuning is larger there. 

We note that the least tuned type I squark models can achieve Apx as low as lO'^, and we stress 
that all of these models have regions of parameter space that are significantly less tuned than all type 
I Higgs and nearly all type II models. 



- 17 - 




0.1 0.2 0.3 0.4 0.5 0.1 0.2 0.3 0.4 0.5 0.6 

A/M A/M 



Figure 3. Examples of the type I squark model parameter space. Left: Qips d4'5 l with Nm ~ 6. Right: 
U'hi,D'H,D with Nm = 3. In the left (right) model, QsiUs) couples to the messengers, so above the horn 
feature, UsiQs) is tachyonic (3.6). To the right of the horn, Q3{U3) is tachyonic (3.12). In every type I squark 
model, the underside of the horn is the region of least tuning. A second region of comparably low tuning sits 
in the middle of the horn for these two particular models. 



3.3 Type II couplings with mixing 

There are five type II (MSSM-MSSM-Messenger) couplings where the messenger mixes with one of 
the fields in an MSSM Yukawa coupling. Three of these couplings are top- Yukawa-like, QU(t)5,H^,, 
UHy^(j)iQ Q and QHu4'w,u ^ ^-nd two are bottom- Yukawa-like, QH^icj)^ and QD(j)c^ fj^. The three top- 
Yukawa-like models are especially interesting because they provide an additional enhancement to At 
by allowing two fields in {iJu, Qs, C/3} to contribute to At, (1-4). This enhancement is so effective that 
one of these models, UHu4>io,q, is the least tuned of all models, possessing regions with Apx ~ 850. 
Overall, all three of the yt mixed models are significantly less tuned than the other type II models, 
because of this At enhancement. 

In order to present the formulas for any type II model (with a coupling of the form XXiX2(l)x-)) in 
a general simplified form, we first define di = d^^"^^ (similarly for ^2, d^). We can derive an expression 
from (2.23) which can be applied to each type II model - mixed or unmixed.'^ For these models, again 

^However, the QQ4ijy can not be treated with these formulas due to the repeated Q. As this model is one of the 
most tuned models, we will not address it in detail. It is straightforward to derive the soft parameters in that model 
from (2.23). 
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# 


Model 


di 


rf2 


rf3 


C7- 


III 




1 


2 


3 


(13 3 8\ 
V 30 ' 2 ' 3 / 


II.2 


UHu4>w,Q 


2 


3 


1 


(13 3 8\ 
V30 ' 2 ' 3/ 


II.3 


QHu4'w,u 


1 


3 


2 


('13 3 8\ 
V 30 ' 2 ' 3 / 


II.4 




1 


2 


3 


(' 7 3 8\ 
V 30 ' 2 ' 3 / 


II.5 


QHd(h,D 


1 


3 


2 


(' 7 3 8\ 
V 30 ' 2 ' 3 / 


II.6 


QQ'^5,D 


2 


2 


4 


f- 5 4) 

VlO' 2'^/ 


II.7 




2 


2 


2 


(i'0,4) 


II.8 




1 


3 


2 


7 3 8^ 
V 30 ' 2 ' 3 / 


II.9 




1 


3 


1 


(-14 8) 
Vl5'"' 3/ 


11.10 


HuD(t)24,X 


3 


2 


1 


3 8\ 
V 30 ' 2 ' 3 / 


11.11 




1 


1 


2 


\ 10 ' 2 ' / 


11.12 


HuL(f>24.S 


1 


1 


2 


V 10 ' 2 ' 


11.13 


HuL4>2i.w 


3 
2 


3 
2 


1 


(to' 2 ' 


11.14 




1 


1 


2 


V 10 ' 2 ' 


11.15 


HuHd4>2i,s 


1 


1 


2 




11.16 


HuHd4>2A.W 


3 
2 


3 
2 


1 





Table 4. The definitions of the parameters appearing in the formulas for the type II models. The first five 
entries are the type II models with MSSM-messenger mixing. The di values preserve the order of the fields in 
the model column, e.g. for model II. 1 di = d^"^"" , d2 — d^**"" and da = d^^- As before, Cr ~ c^" +c$^ +cf^ 
is the sum of the quadratic Casimirs. Note that model II. 6 can not be directly plugged into the formulas of 
(3.13). 



with coupling XXiX2(l)x3, we have, 

Smj,^ = (^d, ^ d,\^ + 2d^d^\^yl2^ - 2d^Cr9lX^ 

- dl^d2X-yl2, + \d,d^2'XV2,,) - '-^^^^2 



(3.13) 



Smj^^ = Smj^^{l o 2} 

Sm'x^ = (dl'd^yl^ + d'j'd2yl,) X'A' 

A - ^1'^ \2 a 
^^--"16^^"^ 

The 2did^X^yf23 piece of Sm'j^_^ vanishes unless there is MSSM-messenger mixing, i.e. both XXiX2(j)x3 
and 2/i23A"iA'2A'3 appear in the superpotential. The di values are tabulated for each type II coupling 
in table 4. We now turn our focus to the individual models. 

QU4>^.Hu — This model is unique in that it is not bounded at high A by any stop tachyons, 
although slepton tachyons provide a similar ceiling. Both stops receive a substantial enhancement 
to their soft masses from the mixing generated y^A^ term. Unsurprisingly, this model occupies a 
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parameter space with features very similar to many type I squark models. However, unlike those 
models, it is least tuned at high A and high just above the region where At is large enough to 
contribute negatively to the m/j. The tuning in this model is controlled by A a and /i at high and low 
jj:, respectively. The tuning in this model is shown in fig. 5. 

UHu(t>w.Q — Due to contributions from both Au and Ah^, this model receives the largest At 
of any type II model. At — j^^- Additionally, this is the least tuned of any model ~ type I or II. 
The tuning contours are shown in fig. 4. This model possesses many aspects of both the type I Higgs 
and squark models. Qa is tachyonic at a relatively low A, however, the U3 and tachyons intersect 
near (A, ^ (1.2, 0.32), so no horn feature appears in this model. As with the type I Higgs models, 
this model is least tuned for small because the one loop contribution to m?^ is large everywhere 
else. This leads a strip in near A ~ 1 of lowest tuning. This strip is cut off from above by the 
large positive contribution to rn?fj^ causing problems with EWSB (as in type I Higgs models), grows 
more tuned to the right from the increasingly negative rii^H^ more tuned below by an increasing 
messenger scale. In the region of lowest tuning, A3 is the largest contributor. For high and high A, 
tuning is controlled by the one-loop contribution. At low A, /z is the largest contributor to tuning. 

QHu4>w,u — This model, shown in fig. 4, is similar to UHu(f>io,Q, however the relation du — 2dQ 
works against this model in two ways. First, the smaller multiplicity factor for Q3 means that the 
A-term is not quite as large as in the UHu4>io,Q model. Second, the larger multiplicity factor for 
gives a large contribution to from the mixing term y^A^ and this contribution causes problems 
with EWSB at lower A compared to the UHu4>io.q model. In particular, this cuts into the region 
where At/ Ms — > a/6, further reducing the quality of this model. As in the UHu4>io.q model, the 
tuning at low A is controlled by /i. At higher A values, A3, A^ and Ai_/oop control the tuning for low, 
medium and high respectively. 

QHii(j>j^ — This model lives in a short slice of parameter space bounded by U3 tachyons above 
and below by regions where EWSB cannot be achieved. The tuning in the model is controlled by 
everywhere. Since the mixing of this model gives additional enhancements proportional to y^, the 
effect of the mixing leads to only a minuscule contribution which makes no appreciable change to the 
case where the term is absent. 

QDcj)^ — This model lives in a rather narrow slice of parameter space bounded by tachyons 
above and D3 tachyons below. The tuning in the model is governed everywhere by fj, and A^. As the 
in the previous model which mixes with the bottom Yukawa, no significant change manifests from the 
effect of the mixing. 

3.4 Type II couplings without mixing 

The type II couplings without mixing, which frequently have very small viable regions in parameter 
space, tend to have high tuning. This is primarily because near a stop tachyon, Ms is smaller (allowing 
for large At/Ms)^ but the viable regions in many type II models are sculpted in part by tachyons which 
are uncorrelated with Ms- Additionally, these models do not receive any significant At enhancement, 
as was the case for models II.1-II.3 as discussed in the previous subsection. We now briefiy discuss 
the features of each model. 

QQ4>b,D — Naively, one might expect this model to significantly enhance At and manifest with 
regions of low tuning. However, it turns out that this model is prevented from ever achieving maximal 
mixing. First, since Q3 contributes doubly to Ti^^, this results in a maximum MSSM-messenger 
coupling of A ~ 0.8 before C/3 tachyons enter. Additionally, Qz itself is tachyonic for all between A ~ 
{0.45, 1.1}. This is due to the same contributions which induce a squark throat (i.e. the discriminant 
in (3.12) is negative). Thus, this model has no valid solutions in any region with A > 0.45. The tuning 
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Figure 4. The tuning shown for the two least tuned type II models, both of which have mixing between the 
messenger and an MSSM superfield. Left: U Hu4>\o,Q- This model has mixing between the MSSM superfield, 
Q, and the messenger 4'io,Q- Right: QHn<j>iQ^u ■ This model has mixing between the MSSM superfield, U , and 
the messenger 4>w,u- In both models, the tuning increases at high -p because of the rising 1-loop contribution 
to TO^^ . Both parameter regions are bounded on the right by stop tachyons. Above, they are also bounded 
by stop tachyons at higher while at lower the issue is unattainable EWSB due to the large positive 
contribution to . 



in the small region that is valid in this model is governed by /i and A3, however, the MSSM-messenger 
coupling contribution is truly negligible here. 

UDcj)^ jj — This model, with tuning shown in fig. 5, lives in a small slice of parameter space 
bounded by Q3 tachyons above and U3/D3 tachyons below. The regions with the best tuning in this 
model are at higher values of A near where is tachyonic. Encroaching on the region bounded below 
drives the SUSY breaking scale A up rather drastically due to the At /Ms growing so large that it 
provides a negative radiative correction to (this effect appears so pronounced in this model mostly 
due to the small parameter space considered). Near the very large A regions, the tuning is set by /i, 
but in the regions of lower tuning, it is set by the A^ term. Of all the unmixed type II models, this 
one presents with the lowest tuning. 

QL(j}-^ ]j — Perhaps unsurprisingly, this model is very similar to the QHdiib model in the previous 
subsection. It lives in a short slice of parameter space bounded by tachyons above and by slepton 
tachyons below. Here, the tuning is controlled by A^ and /i. 

U E(f>^ — This model is narrow in A, but stretches further into ^ than the other type II models 
before the Q3 tachyons above, and the f/j tachyons below close the region. The tuning in the model 
is controlled by fi everywhere. 

HuD(t)24^x — This coupling involves ff„, D and a superfield with gauge charges like the X bosons 
of Pati-Salam models. This model has similarities with both the Higgs and squark type I models. At 
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Figure 5. The tuning shown for the third and fourth least tuned type 11 models,. Left: QU (fi-g . This model 
is the worst of the top- Yukawa-like mixed models. Unsurprisingly, it also receives the smallest contribution to 
At- It generally looks similar to the type I squark models in shape. The tuning is best however at high 
Right: UDcj>c, jy. This model is bounded from above by tachyonic stops from Q3. The divergent region at the 
bottom of the plot is where At/ Ms is growing non-perturbative. While much of the region has lower tuning, 
it exists in a very small window of A and space. 



low , raising A causes issues with electroweak symmetry breaking. To the right a slight tachyon 
throat cuts off the model. The tuning is driven by /i, Ki^ioop and A^. 

HuHd4>i,s^ HuHd4'2i,s, HuHd4>2i.w — In addition to being horribly tuned, these models introduce 
a ^ — Bji problem. Overall, the other features of these models is quite similar to the type I Higgs 
models. Note that for H.aHd(f>i.s, the singlet field is insufhcient to mediate SUSY breaking, so an 
additional 5-1-5 with no superpotential interactions with the MSSM is assumed. The tuning is driven 
by Ai_;oop and fi. 

HuL(j)i^S, HuL(j)2i^Si HuL(l)24.w — These models have EWSB issues above and a slepton tachyon 
throat structure analogous to the squark throat discussed in section 3.2. However, more importantly 
these models introduce unacceptably large neutrino masses. Additionally, the (jji singlet field is insuf- 
ficient to mediate SUSY breaking, so an additional 5 -f 5 with no superpotential interactions with the 
MSSM is assumed. The tuning is driven by Ai_ioop and /i. 

4 Phenomenology 

In this section, we discuss aspects of the phenomenology of the models that give Aj^j- < 2 x 10'^. The 
models that satisfy this are: all type I squark models, the three mixed type II models - U^Hu4>io,q, 
Q'iHu(t>io,u and Q3U3(p5^H^ ~ and one other type II model - U^D^cf)!^ jj. Nearly all of these models 
possess spectra just beyond the reach of existing LHC searches. It is interesting to note that such heavy 
spectra actually seem to be a requirement of these GMSB models with a 125 GeV Higgs. This strongly 



- 22 - 



1.8 1.9 I.IO I.ll 1.12 1.13 1.14 1.15 1.9' 1.13' II.l II.2 II.3 II.7 
5000 ^ ^ ^ 5000 



3000 



^ 1000 
p. 700 



1 1 1 1 1 1 1 1 1 1 













- ^— ~ — — — ~ 


— 


~ ~0 ~ + 

^ X X 

^ 




1 1 1 1 1 1 1 1 1 1 


\ \ \ \ 



250 



100 

1.8 1.9 1.10 I.ll 1.12 1.13 1.14 1.15 1.9' 1.13' II.l II.2 II.3 II.7 

Model 



3000 
2000 

1000 

700 

500 

250 
100 



Figure 6. The spectra for some of the better models at their points of least tuning are shown. All type I 
squark models are shown to the left (Q: 1.8-11 and U: 1.12-15), type II models, including the three models 
which mix the top Yukawa with the messenger field and the UD(j)D (II-7) are shown to the right. 1.9' and 1.13' 
denote the best point within the distinct region of comparable tuning accessible in these two models (see fig. 3) 
which present a very different spectra. In the plot, thick, large lines denote colored particles - g, t\, t2, 6i, 62 
and q (the nearly degenerate first-generation squarks) are shown. The thinner lines denote uncolored particles 
- I, x° Eind X* s-re shown. All four neutralinos and both charginos are displayed. In nearly all models, all 
right-handed sleptons and all left-handed sleptons/sneutrinos are approximately degenerate. 



suggests that the non-observation of SUSY and the presence of a heavy Higgs may be correlated issues 
rather than two distinct problems of SUSY. 

4.1 Type I squark models 

In the region of least tuning (the base of the horn in fig. 3), the type I squark models have heavy 
gluinos and first generation squarks falling between ^3.5-5 and ~3-4.5 TeV respectively, while the 
lightest stop (as well as the sbottom in models) has a mass between ^0.5-1 TeV. Additionally, 
there is almost always an NLSP f or co-NLSP ^s generally between ^300-500 GeV (although these 
sometimes appear even heavier than 700 GeV). However, the other region of low tuning appearing in 
models L9 and 1.13 (in the center of the horn) has a rather different profile (the best points of this 
second region are denoted by L9' and 1.13' in fig. 6). Here, the models have heavier stops, ^1.2-2 
TeV, but since A has dropped significantly, the gluinos and first-generation squarks are now much 
lighter ~2.0-3.5 TeV and ~1.5-3 TeV, respectively. Surveying these points with less tuning, it is clear 
that the mass of the lightest stop and the masses of the gluino and first-generation squarks tend to be 
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anti-correlated ~ thus, much of the parameter space with lower tuning possesses either squarks/gluinos 
or stops which can accessible at 14 TeV LHC. 

Given a stop portal, there are three separate simplified topologies that arise: 

1. t "NLSP" — with the stop appearing as either a co-NLSP with £fi, or as the true NLSP (e.g. 
the best point of model 1. 11 in fig. 6) 

2. i : B : e — Decays of i tx" ^ U'^i"^ 

3. i : i — This has competing decays of f — > biyf^ through an off-shell H and i — > t£^l^ 

In all three cases, the NLSP decays to its SM partner plus gravitino. The first case is rare in these 
models, but has been discussed in the literature extensively, see e.g. [34-38]. The second case, which 
is slightly more common, will populate high I^t, high Ht multilepton searches ~ bounds for these 
will likely fall at or near the kinematic limit for titl production. The third case happens in models 
1.8-10,12-15 in fig. 6. The third case is most interesting in the limit where the top decay is squeezed 
out (toj — TTif < rrit). A sample spectrum which produces this signature is given by the best point of 
model 1:14, as shown in fig. 6. 

In this third case, when only the Higgsino-mediated decay can occur, the signature is two 5-jets 
with moderate px, opposite sign rs with very high px and large additional I^t- While in principle, 
several existing searches could have sensitivity to this exotic signature [39-42], those searches are 
not optimized to handle this kinematic configuration. Existing opposite sign r -I- I^t searches tend 
to require a very hard jet and/or high Ht- In the squeezed regime, the 6-jets will generally lack 
substantial enough pT to meet these strict cuts. Third generation leptoquark searches would have 
some sensitivity, however, the br invariant mass and St requirements may be too harsh for this 
specific signature since so much of the event's energy is put into I^t- Searches for Mt2 in the tt 
system are likely the most sensitive probe of this system. Estimation based on the cuts in [42] (as 
done for the similar {hr^ vv}{bT~ vv} signal examined in [43]) suggest this study repurposed would be 
sensitive to stops in the 450-550 GeV range. A dedicated search region either for rs in the Mt2 search 
or a hard pt,t region with softer cuts on the jets could likely probe significantly higher stop masses. 
Further exploration of this signature is an exciting subject for future work. 

The distinct regions of low tuning appearing with heavier stops and lower A (in models 1.9 and 
1.13) can have accessible gluinos and first-generation squarks (these spectra are denoted 1.9' and 1.13' in 
fig. 6). These production portals result in complicated SUSY decay topologies with slepton co-NLSPs. 
These models will generically give hard jets, large I^t, and multiple leptons. The 14 TeV extension 
of existing SUSY searches will be very sensitive to these topologies if the colored sparticles are light 
enough for a meaningful production cross-section. Regions in other type I squark models exhibit the 
same decreasing A, but pay a substantial price in tuning. While the spectra are not shown, they 
manifest with the same qualitative features - rising stop masses with all other sparticles decreasing - 
which allow for first-generation squark/gluino production. 

4.2 Type II models 

As one can easily infer from fig. 6, the type II models have more variation in their signatures. Generally 
these models possess lighter gluinos and first generation squarks. In fact, the points of best tuning 
in model II. 1, the mixed QsU^c/j^^h^ model, are already ruled out by existing searches because these 
colored particles appear near or even below 1 TeV. However, in the other three less tuned models of 
type II, the gluinos are closer to 2 TeV and are not necessarily constrained by existing searches. 
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The least tuned type II model, II. 2 {U^Hu(f)w.Q), has stops and bottoms near 1.5 TeV, but the 
other squarks and gluinos are near 2 TeV. For very low these models will often have a B NLSP, 
however, the bulk of parameter space does have a right-handed slepton co-NLSPs. In these low jj 
regions, the SUSY breaking scale can be high enough to for a detector stable neutralino. These models 
can present with classic SUSY signatures of multiple leptons, jets, Ifx and very high St- 

The colored states of model II. 3, QsHucfiio^u, are slightly heavy, with stops near 1 TeV and both 
first-generation squarks and gluinos appearing near 2.5 TeV. This spectra should be observable at 
14 TeV LHC, where the production cross-section for 2.5 TeV gluinos and squarks is O (l fb~^) [44]. 
However, even moving into regions with slightly higher masses, these models still have very light 
slepton co-NLSPs appearing near 200-300 GeV, which could be easily observed at a future ILC. 

A slightly different topology, with stops, gluinos and first-generation squarks near 750 GeV, 1.5 
TeV and 2 TeV, respectively manifests in model II. 7, U^D^cf)-^ jj. Here, the Higgsino is very heavy, 
but the O (250 GeV) bino is always the NLSP. As above, these complicated topologies cannot be seen 
until 14 TeV, but give rise to classic SUSY signatures of multiple leptons, jets, Ifx and very high St- 

5 Conclusions 

In this work, we studied models that produce large A-terms through the introduction of a single 
marginal superpotential interaction between MSSM fields and the messengers of minimal GMSB. We 
classified all such interactions compatible with perturbative SU{b) unification. Our complete list of 
31 possible couplings - 15 type I couphngs (MSSM-messenger-messenger) and 16 type II couplings 
(MSSM-MSSM-messenger) ~ is summarized in table 1. 

Motivated by rampant confusion in the literature concerning the correct soft term contributions 
from MSSM-messenger interactions in the presence of MSSM-messenger mixing, we derived a new 
method for treating these interactions by directly integrating the wave functions in a manifestly holo- 
morphic scheme. This conceptually straightforward method produced results applicable to all scenar- 
ios, whether or not fields are mixed. Formulas were presented for the soft parameters both in complete 
generality (2.15); and in the simpler special cases of only type I and type II couplings, (2.19)-(2.20) 
and (2.22)-(2.23) respectively. 

Using these new formulas and a slight variation on the Barbieri-Giudice tuning measure, we sur- 
veyed the tuning in each of the 31 models. Under this examination, we concluded that the qualitatively 
similar type I Higgs models universally have high tuning due to the little A/niH problem, while the 
type I squark models can have tuning as low as ApT ~ 10'^. The majority of type II models have 
poor tuning, with the notable exception of the three models which allow for MSSM-messenger mixing 
with a top- Yukawa-like interaction, which generate a very large value for At- The least tuned of these 
models, U^H^ipQi manifests with the lowest tuning of any model studied in this work, /S.ft ^ 850, 
while the other two still have tuning below Apt < 2 x 10'^. 

The spectra in the least tuned models usually have particles beyond the reach of the 8 TeV LHC, 
but can be accessible to the 14 TeV upgrade. This was not put in by hand, but is a consequence 
of the requirement that rn^ = 125 GeV, together with the minimal gauge mediated structure of the 
messenger sector. This suggests that the failure to find superpartners so far at the LHC was not an 
accident, but in fact had to be the case. 

The models tend to possess either a stop below a TeV (and sometimes a sbottom as well) or 
accessible first-generation squarks and gluinos. The decay chain in these scenarios terminates with 
the light NLSP stau, co-NLSP sleptons or NLSP bino (with the rare exception of a stop NLSP) 
decaying to a gravitino. Both prompt and long-lived NLSP decays are possible in the models that 
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we have considered. Production of the first-generation squarks and giuinos gives rise to classic SUSY 
signatures of high I^t and St with muhiple leptons and hard jets. The scenarios with only stop 
production will usually given multi-lepton signatures; however, the cases with a stop NLSP or a 
squeezed t ^ f^bv transition are more difficult. The latter gives rise to a very interesting and poorly 
studied bbT^T~ + I$t signature. Uncovering search strategies to improve sensitivity to this final state 
is an exciting avenue for future study. Another very common feature in these models is light sleptons 
(and occasionally Higginos) which could be readily produced and studied exhaustively at a TeV scale 
ILC. 

While we took an agnostic approach to flavor physics in this work (assuming a perfect alignment), 
the lack of MFV in all of the least tuned models begs a thorough treatment of flavor. How much 
misalignment is permissible and whether alignment can be naturally achieved in a sensible way are 
both questions for future study. Another issue that we have not addressed in this paper is the origin of 
/i and -Bp. It would be very interesting to extend the models considered here to include a mechanism 
for n/Bfj^. Perhaps extensions involving the NMSSM along the lines of [21] are viable. Finally, while we 
assumed single MSSM-messenger interaction terms for simplicity in this paper, it would be interesting 
to explore the effects of including of multiple MSSM-messenger interactions, as this could potentially 
generate regions of parameter space which exhibit even lower tuning. 
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A A Detailed Study of the QU^ Model 
A.l Applying our general formulas 

In this appendix, we will provide an in-depth study of the model with 

W = XQU<^ + ytQUH^+X<i>^ (A.l) 

with {X) = M + 6^F. (We drop the 3rd generation subscript on Q and U to avoid cluttering the 
equations.) In this model, the messenger $ has the same quantum numbers as Hu- This is in many 
ways the prime example of a mixed type II model, given that it has been studied already in many 
papers [13-16, 18-20]. We will use this example to illustrate a number of points. First, it will highlight 
some interesting features of our general formulas. Second, by computing the soft masses in this model 
using other methods, it will provide a detailed check of our general formulas. Finally, we will use this 
example to illustrate the shortcomings of the formulas and approach in [12] which were mentioned in 
the body of the paper. 

As seen in (2.23) and (3.13), the effect of MSSM-messenger mixing appears in the |A|^ \yt\^ terms, 
so to focus on that, let us for simplicity set the gauge couplings and all other MSSM Yukawa couplings 
to zero. We begin by quoting the result of the general formula (3.13) for this model. Taking Xi — Q, 
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X2 = U , — Hu, and setting Qr — Ub — Vt — 0, we have: 



Au^--^du\\\^K 

1 



5ml - ——^[dQ{du +dQ + d^m" + 2dQd^\XWt?)K^ (A.2) 
5ml = .^{duidu + dg + d^)\X\^ + 2dud^\X\'\yt\^)A' 
5m\^ = -^M^Q + dv)\X\^\ytft^ 

with no contribution to other soft masses. Here, as in section 3, dq is shorthand for dg*, etc. In the 
MSSM, we have dg = 1, du = 2 and c?$ — 3, but it will be useful to leave these multiplicities general. 
Notice that the last two \X\^ \yt\^ terms from (3.13) have cancelled out of niQ and mfj, leaving only 
the last line induced by the MSSM-messenger mixing. One can also check that by substituting the 
standard beta functions and anomalous dimensions into the formulas of [12], one misses these extra 
terms. In the following subsections, we will study these extra terms in more detail. We will compute 
the soft mass-squareds in this model in two different ways: directly using SUSY correlators as in [22], 
and directly using wavefunction renormalization in the interaction basis. The latter will also illustrate 
the subtleties of wavefunction renormalization which the method derived in section 2 avoids. 

A.2 Calculation using SUSY correlators 

Let us calculate the m| directly, using a supersymmetric correlator formalism along the lines of [22]. 
As in that paper, we separate out the A-term-squared contribution to mf coming from integrating out 
the auxiliary field Fi, and we will denote the remainder by the hatted quantity 



2 - 2 

m,, — rrii 



(A.3) 



This formulation is much simpler computationally, because it allows us to avoid various subtleties 
resulting from the treatment of contact terms and total derivatives. 
The two- loop ©([i^l^) contribution to ttiq is given by 

ml = j d'x2 . . . d\e (Q'(C/(A<i> + yti/.))iQ'(C/^(A*$t + y^Hl))^ 

Q2(Q[/(A$ + ytHu)hQ\Q^U\X*<^^ + y;Hl)U<^^U<f^^^)e) 

= -\F\'dMdM' J d^X2 . . . d^xi (((/(A* + 2/*i70)i(C/t(A*$t + y*Hi))2 (^-4) 
Q^(QC/(A$ + ytH^jhQ\Q^U\X*'f^ + ylHl)),) 

= ~\F\^dMdM'Zg\M,M*) 

These correlators are evaluated in the Euclidean supersymmetric free theory and only contain IPI 
diagrams with respect to the theory containing the auxiliary fields Fq^u^h^ (in other words, if a diagram 
would become disconnected were an auxiliary field propagator removed, then it is not considered IPI.) 
The subscripts i = 1,2,... are shorthand for the positions Xi. In the second equation, we have rotated 
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the supercharges so that they act on the operators at 2:5 and xg, and then transformed this into 9m 5m* 
of a simpler correlator. Comparing with (2.4), we can see that the correlator being differentiated is 
the 2-loop contribution to the wavefunction of Q. This is also clear diagrammatically, if we view this 
as the two-loop IPI diagrams with external FqFq in the presence of the interactions (A.l). 

It is straightforward to expand out all the terms in Zq and perform the free-field contractions 
(keeping in mind the IPI condition). The result is: 

= -\\\'dQ{duh + d^h) - \\\''\yt\''dQ{du{h + h) + 2d^h) (A.5) 



where 



I p^{p^ + \M\^)qmp + q)^ + \M\^)' j {p^ + \M\^yq^{p + qY 



(A.6) 



pV((p + ?)2 + |M|2)' " J p-^{p-^ + \M\^)q^{p + qY 

(27r)4 (27r)4 



with J shorthand for the integral over Euclidean phase space / (2J)4 • easy to check that 



We see that and 1^ contribute with equal magnitude and opposite sign to dMdM" Zq. Thus the 
|Ap|j/(P contribution proportional to dfj drops out, consistent with what we found in (A. 2) using 
the general formula. However, the contribution from mixing proportional to d^ remains. Combining 
(A. 4), (A.5) and (A. 7), and adding in Aq — -jj^^iyi j^fp , we find perfect agreement with (A. 2). 

Next, consider the soft mass for U. This can be found by the same manipulations used to derive 
rriQ, but with dg ^ du- This agrees with the formula for mfj in (A. 2). 

Finally, we come to mjj . Here there is no A-term-squared contribution, and we have: 

m%^^ = -\F\^ J d^X2 . . . d^X, (Q2(2/*Q?7)iQ'(2/*Qt[/t)2g2(^gf;^)^g2(^*gt[/t<i,t)^($$)5(<i,t$t)g 

= -|A|2|y,|2|F|2aM9M* / d^X2 . ..d^x^ {q\QU){q\q^ U^)2{QU <^UQ^ <^>%) 

J \ / (A.8) 

= +d^{dQ + du)\X?\yt?\F\'dMdM'h 

d^idQ+du)\^yUF\' 

(167r2) |M|2 

which again agrees perfectly with (A. 2). 

A. 3 Analytic continuation method 

As another check of the general formulas, let us also perform the analytic continuation calculation by 
directly and explicitly integrating the wavef unctions, which can be easily done in this simple example. 
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The trick is to do a unitary field redefinition to go to the interaction basis: 



V + yt \/A2 + yf 

V + yt V ^ + yt 



So we will study the equivalent theory defined at the scale A: 



= AQL/$i + + K2X$2$ (A.ll) 

K = Q^'Q + U^U + ^{<i>i+<^>l<p2 + ^^^ (A.12) 



where 

A , yt 



A = -^A2 + y2 ^ Ki = ^ , and K2 = (A.13) 
A A 

In the interaction basis, <I>i will receive wavefunction renormalization, but $2 will not. This fact 
simplifies the calculation considerably. 

Now we evolve this theory down to a lower scale. As in the previous subsections, we again neglect 
gauge couplings and the other Yukawas for simplicity. The result is: 

W = XQU^i + kiA:$i$ + K2X^2^ 

K = ZQ{t)Q^Q + Zuit)U^U + + $|$2 + ^^'"^^^ 

where the wavefunctions obey the beta ftmctions: 

^- 1/2 1/2 1/2 

for i = Q,U,^i. Here A(t) = XZq Zjj Z^ is the nonholomorphic Yukawa coupling in the 
theory with canonical Kahler potential. It obeys the beta function 

^=A(t)(7Q+7a + 7*) (A.16) 
with boundary condition A(0) = A. These equations can be integrated to obtain Zi{t)\ the result is 

= dq + du + d, 1^^ 8^^^ ') ^^-^'^ 

Next, at the scale M, we integrate out the messengers supersymmetrically. This sets $ = and 
$2 = —1^1^1/1^2- We additionally identify $1 with Hy, below the messenger scale. At M this theory 
then becomes 

W = XQUHu (A. 18) 

K = ZQ{tM)Q^Q + Zu{tM)U^U + (^Z^,{tM) + ^jHlHu (A.19) 

where Im = log \M\/A. It is clear that the wavefunction of is discontinuous at the messenger scale. 
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This is precisely the issue alluded to in section 2 arising from the mixing of and $ that prevents 
the formulas of [12] from being applied as presented. 

Finally, we should proceed with wavefunction renormalization by again integrating (A.15)-(A.16), 
but now from M to some lower scale fi. The key difference is that because of the discontinuity in the 
wavefunction, the boundary condition for the Yukawa coupling is also discontinuous: 



Taking this into account, we find 



(A.20) 



log^Q,t/(logM) 
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-logf 1 



log 1 
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, 87r2 
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log 
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log Zh^ (log /i) = - log 1 - 



du + rf<E,)A2 , \M\ 



(A.21) 
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log 
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Differentiating these expressions with respect to M and M* and substituting (A. 13) to recover the 
dependence on the original couplings, we again find perfect agreement with the general results (A. 2) 
for w^,t/,ff„- 



B Fine- Tuning Measure 

Fine-tuning is an inherently ambiguous concept. When comparing variations of two nearly identical 
parameters, it seems sensible, however even when comparing the variation of a mass parameter to a 
coupling the choice of measure quickly looks somewhat arbitrary. Since one of the objectives of this 
work is to determine which GMSB models possess lower fine-tuning, we demand that our fine-tuning 
measure, A.ft, satisfy certain properties: 

1. ApT should provide a meaningful and accurate comparison between GMSB scenarios 

2. ApT should never overlook contributions from large terms which cancel in a uncorrelated way 

3. ApT should never introduce contributions from large terms which cancel in a correlated way 

4. ApT should assign comparable sensitivity to two uncorrelated terms which cancel one another 
Traditionally, the Barbieri-Guidice tuning measure [26] is defined as: 

r A 1 1 A d log m'j /-r. , N 

Abg = niax{Aa} where = ° B.l 

a log a 

where a sums over a set of "fundamental" parameters. In general, it is not so clear what these 
fundamental parameters should be, and different choices for them lead to different numerical values of 
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the tuning measure (the measure is not reparametrization invariant). For our purposes, we will adapt 
the Barbieri-Giudice tuning measure to be, 

d log Tn?' 

Aft = max{A,} where A, = " (B.2) 

o log A| 

where Aj G {A^A, Ai_/oopj /^} and runs over all the important couplings in the theory, i.e. in 
this case, A^ £ {51, 52, 53, J/t, J/6, J/r, A}, although in practice only variations in (73, yt or A mani- 
fest deviations large enough to matter quantitatively for this study. Thus, our set of parameters 
is Ai € {A3, At, Aa, Ai_ioop, /i}. As in [45], we choose to differentiate ml in (B.2) with respect to only 
parameters with mass squared units as this serves to better adhere to our requirement that canceling 
terms provide comparable sensitivities. 

With the exception of /i (which is calculated directly), we compute these derivatives by imple- 
menting a very small fractional change in the parameters injected at the messenger scale, run down 
to the low scale and measure the change in ml using SOFTSUSY. Ai^ioop, see (2.16), was chosen as 
the parameter to account for the dependence on (i.e. we fractionally vary the term j^\J h (-^) 
rather than varying M, or some other combination of these parameters.). To keep variations 
in Ai_;oop and A^ orthogonal, we keep the 1-loop term fixed when we vary A^. This method removes 
the possibility for uncorrelated cancellations and correlated values being erroneously treated. Addi- 
tionally, it is well defined in all of the GMSB scenarios we present here and, in principle, can be easily 
translated to work for many other GMSB scenarios as well. 
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